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1. Introduction 

Nanoelectromechanical systems (NEMS) are systems where mechanical and electronic 
degrees of freedom are coupled and whose characteristic length scales are measured 
in nanometers. While metal or silicon are common material choices for 
microelectromechanical devices (MEMS), carbon nanotubes (CNT) may become one 
of the mainstays in future NEMS technology [1] due to their unique combination of 
electrical and mechanical properties: small mass, extraordinary stiffness, low mechanical 
dissipation and electrical properties ranging from semiconducting to conducting [2]. 
Combined together, these properties allow for NEMS devices operating in the high GHz 
regime. Several CNT-based NEMS devices have already been demonstrated [31 HI [5] . 

In [HI [7J [8] one specific such carbon-nanotube-based system was considered: a three- 
terminal nanomechanical relay with a layout similar the one in figure [U Such devices 
have since been successfully fabricated [H [101 [11]. In [6] [9] [10], E] the device was 
considered mainly operating as a switch and transduction was based on a tunneling 
current between the tube tip and the drain terminal. In a subsequent publication [12] 
also the high frequency properties were considered and nonlinear resonant behavior was 
demonstrated. 

Nonlinear response to AC-drive is a characteristic feature of many MEMS devices 
as well as in NEMS p3] and has recently received much attention [Til [T5l fT6[ fT7 [ fT8 l fT9] . 
It can typically be modeled by using the Duffing equation. Understanding of the 
parametric dependence on the behavior of this response is essential for any potential 
technological application. Of equal importance for technological applications is to 
understand how NEMS devices act when embedded in an electronic circuit. In this paper 
we investigate more carefully the RF-response of the three-terminal CNT resonator 
in figure [TJ subjected to a harmonic input signal on the gate electrode, going beyond 
the simple Duffing equation. In contrast to previous publications [32], we base the 
transduction not on electron tunneling between the tube tip and the drain but on the 
displacement current generated by the time-varying tube-drain capacitance when the 
device is connected to lossless transmission lines. 

This paper consists of several sections. In section 2 we present the system and 
derive a lumped dynamical model. Then, by means of perturbation theory, we derive in 
section 3 a set of algebraic equations that can be efficiently used to obtain and classify 
the frequency response of the device along with expressions for the output power (S21- 
parameter). In section 4 we consider the regime of linear response solutions and in 
section 5 we characterize the nonlinear behavior. Finally in section 6 we discuss the 
domain of validity of the perturbative approach by direct comparison with numerical 
integration and discuss the delivered power arising from a purely capacitive transduction. 
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Figure 1. Physical layout of a three-terminal carbon nanotube resonator. A carbon 
nanotube or carbon fiber with diameter D is mounted on a support - the source contact 
- of height H extruding a distance L out over the substrate. Beneath the tube are two 
metallic strips - the Gate and the Drain - of widths Wg and Wd respectively. The 
distances between the gate and drain strips from the support are Zq and Zd and the 
strips have heights Hq and Hd respectively. 



2. Model 

The system is depicted in figure [TJ Mounted on a conducting support (Source) of height 
if is a nanotube of length L with outer diameter D Q and an inner diameter D\. Below 
the tube two conducting strips of heights Hg,d an d widths Wq,d act as gate and drain 
electrodes. The motion of the carbon nanotube can be described as a simple elastic 
beam deflecting in only one direction [20] using the action 

Sbeam = t dt ^ dx^[d t u(x,t)] 2 - Sc^X, t)] 2 . 
Jt x JO I I 

Here u(x,t) is the instantaneous deviation of the tube towards the Drain electrode, at 
a distance x from the tube support. E is the effective Young modulus [201 12T1 122] of 
the beam and I = 7t(-Dq — Df)/64 the moment of inertia. The cross-section area is 
A = \n(Dl — Df) and p is the density of the tube. The external forces acting on the 
tube, i.e. actuation and transduction, arise from capacitive coupling between the tube 
and the electrodes. These forces depend on the instantaneous charge distribution and 
geometrical configuration of the tube. This gives rise to an additional part of the action 
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Figure 2. Equivalent circuit model. The gate-tube capacitance Cq as well as the 
drain-tube capacitance Co depend on the instantaneous deflection of the cantilever. 
The gate and drain contacts are assumed to be connected to transmission lines with 
characteristic impedance Zq. 



for the tube which can be written 



f -L \ 

Sel.— mech. = / dt~ / J 



where CcdKs, t), t] are the tube-gate and tube-drain capacitances at time t and Q' G D 
the associated capacitor charges. Stationarity of the action with respect to u(x,t) then 
provides the equation of motion for the beam 

pA$u(x, t) + EId*u(x, f) = £ £ [A^(t)] 2 ^M (1) 

where we have introduced the shorthand 
5Ci[u,t) _ f 5d[u,t] 



Su(x) J 5u(x,t') 
to emphasize that the functional derivative only affects the spatial dependence of the 
capacitance at time t and the potential differences AVg,d = Vt — Vg,d between the tube 
and the gate/drain electrodes. 

The electronic degrees of freedom are most conveniently treated within the circuit 
model in figure [2j For finite source-tube resistance i?g we have three relevant degrees 
of freedom Q G (t) = Q' G (t) + Q GD (t), Q D (t) = Q' D (t) + Q GD (t) and Q T (t) = Q' D (t) + 
Qx(t) — Q' G (t) ~ Qs(t). With only one incoming signal on the gate, V + (t), these charges 
obey the equations of motion 

(2) 
(3) 
(4) 

where we have allowed for constant DC-bias offsets on each of the electrodes and assumed 
the gate and drain connected to lossless transmission lines with a real impedance Z . 



d t Q G = 


1 
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2V + (t)-V G (t) + V£ 


d t Qv = 
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y D (t) - v° 




d t Q T = 


~B-[Vs ~ V T (t)], 
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For the tube diameters and signal frequencies we consider, the kinetic inductance of the 
CNT [23| [21] is insignificant in relation to other impedances. The capacitances give us 
a linear relation ship between voltages (Vq, Vq, Vt) and charges (Qg, Qd, Qt), 

Qg = (C G + C GD )V G - C G V T - C GD V D (5) 
Qd = C B V T + C GD V G - (C D + C GD ) V D (6) 
Qt = (C T + C s + C D + C G )V T - CgVg - C D F D - C G V G (7) 

ensuring that we have a closed set of equations for the dynamics consisting of 
equations (pD)-^. 

The full PDE in ([1]) is not very tractable as it stands. We will assume that only 
the lowest lying fundamental vibration mode is excited by the incoming signal. In this 
approximation the deformation of the tube can be parameterized by the displacement 
ix of the cantilever tip from the static equilibrium position (for details see Appendix 
A) 

x T (t) + Tx T + n 2 x T (t) = -±— m(t)] 2 ci(x T ). (8) 

where M eff = pAL/5.684 and fi = 3.516L~ 2 ^EI/pA. Also, in © a factor Tx T has 
been incorporated in a phenomenological way to account for dissipation [22j 125] . 



3. Perturbation theory 



Although direct numerical integration of the lumped model is a straightforward task 
it is time consuming due to the multiple time scales involved. It is instead our aim 
in this paper to derive a set of algebraic equations which can be used to classify or to 
quickly determine the response of a given geometric or biasing configuration. We do this 
through a perturbative analysis by means of the averaging method [26]. We will thus 
assume that we have a single incoming signal on the Gate electrode V + (t) = V + cos(fit) 
and cast the system into a dimensionless form by writing t — t/Q , xt = £,H, Ci = qCo, 
V% = ViV , Qi = qiV C , 7 = T/Q and to = Q/Q ] 



£ + 7fr + f 



<?g(7J 
9d(t) 



= |[^G(r)-^(r)]V G (0 + |K(r) 
2v + (t) cos cut — vq(t) + v[ 



^(r)] 2 cb(0 



D 



[v s - v t (t)] . 



-1 Zp 

Here we have defined a = VqCq/ 'M c ^VLqH 2 = £c/£k and e 

Co = 



VLqZqCq. For Co we chose 

C G C^ where C G D are the gate-tube and drain-tube capacitances when the 
tube is undeflected. Typically the capacitances in nanoscale devices lie in the attofarad 
range which implies that for a single device the impedance mismatch is large. Hence, 
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f < 1 providing a good starting point for perturbation theory. The linear relationships 
between charges and voltages imply that we can expand in e 

oo oo 

Qi = E^\r), v l = J2e n vl n \r),t = G,D,T. (9) 

n=0 n=0 

To zeroth order in e we then have 

Vq\t) = 2v + COS UJT + Vq 

v ( »\t)=vI v?(r)=v s . 

Inserting these zeroth order solutions into the dynamic equation for the tip motion we 
find 

| + 7 Ct + $ = f [2«+ cos ur - (v s - v° G )] 2 c> G (0 + ^[v° D - v s ] 2 c> D (0- 
We now make the Ansatz of an oscillatory solution with slowly changing parameters, 

£ = x (t) + r(r) cos[cjr + 0(r)], 
and assume that the quantities f/wr, 0/cj0 <C 1 keeping only the lowest order terms 

7^0 + x o — (2ru; + jru) sin(cjr + 0) 

+ r(l — c<j 2 — 2cj0) cos(cjr + 0) = — i^(x , r, 0, r). 

2 

Here the kernel K(x , r, 0, r) is defined as 

K = [2t> + cos cur — (f S — 'Wq)] 2 c / g (xo + r cos[u;r + 0]) 

+ ( v d ~ ^s) 2 c D (a;o + r cos[cur + 0]). 

Provided the slow variables r, and x do not change appreciably during one period we 
can average over one period to find 



7i; + ^o = i—/ dr K(x , r, 0, r) 

47T JO 

— / dr K sm wr + 
4tt Jo 



— -'jruj 



-r(l — cj 2 ) — ro;0 = — / dr fTcos(W + 0). 

2 47T Jo 

For generic capacitances we would not expect to be able to do these integrals exactly. 
However, for the case when + wn/2 < I the derivative of the capacitances can 
be approximated by an inverse square law according to the local tube deflection, 
c 'g,d(0 = Cg, d «g,d(«g,d - 0~ 2 wi th 

a G , D = (1 - H G>B /H)u (L)/u (Z GtD /L + W G , /2L). 

Introducing 



og,d = ^ (cxg,d - x + yj (a G ,D ~ x ) 2 - r^j 
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and v sg = — V G , v sd = — an d "^sg = + "^sg; the results after performing the 
integrals are 



7x0 + xo = 2 ^ c G ttG ^G ^~2 g ( a 2 + -q _ g v+VsGaG cos ^ + 2v 2 + (3 - a G 2 ) cos 20] 

7~ I d(~i J. ) 



, 2<7cg ) a D qg ) (qg ) + 1) <<2 , . 

+ r 2 (4 - 1)3 UsD (1U) 

2fu + 77*0; = 8a „ r f +f sg sin — t> , a c sin 20 (11) 

r («g - !) L ' ' J 

r(l - w ) - 2rw0 = 8a^-^ — [w SG a G - v+v SG {a G + 4 - a G ) cos0 

r (a G — L) 

+v 2 + (a G + 2a G l - a G 3 ) cos 20] + 8a ^\ 3 ^d (12) 

The output signal is the displacement current in the drain contact and we will 
express this in terms of the square modulus of the S-parameter IS21I 2 . In general the 
S-parameter Sij is a complex quantity relating the amplitude and phase of an incoming 
signal Vi on port i to an outgoing signal on port j through = Vi/Vj. In terms of 
power we have then IS21I 2 = Pont I Pin relating the total outgoing RF-power on the drain 
to the incoming RF-power on the gate. For the given input signal v + cos(ujt) the average 
delivered power to the device is P in = Vqv\/2Zq and the total power delivered on the 
output is 

- 

out 2Z n Jo 

From the linear relation between charges and voltages we have 



dr(v° D -v D (r)) 2 . 



qr>{r) = c D v T + c GD v G - (c D + c GD )v D . 
Recalling the perturbation expansion we get to lowest order in e 

v d = v □ + e((v s - ^d) c d(0£ - 2uw + c G d sin^r) + 0(e 2 ). 
Evaluating the integral in the stationary state (xq = = r = 0) one finds 



s 21 \ 2 = A(z c ny 



4 



r 2 (a 2 3 — 1) E 



+ 4 (— )ZtJ2 — COS( f ) + C GB 



v + >r{al - 1) — • ™j ■ (13) 

A similar result can be derived for the reflection coefficient \Sn | 2 - Note that the prefactor 
of e 2 = (Z CqQ) 2 indicates that for mismatched systems (e <C 1) only a very small 
amount of the incoming power is actually delivered to the device and that most is 
reflected back (|5ii| 2 ~ 1). 

In a typical setup we have a G 3> «d- I n this case we have found that the system 
(jl0p - (|12p can be simplified considerably by omitting all terms related to the parametric 



Nonlinear resonance in a three-terminal carbon nanotube resonator 



8 



driving and double frequency components and it suffices to solve the simplified set of 
equations 



■yru 

d- 



aC G ~2 



2a 



G 



J SG 



+ 



2o-c^a; D a^(a^ + 1) 



r 2 (4 - l) 3 



ac, 



G 



«G 



v + vsg sm . 



(14) 
(15) 



CO 



ac r 9 

2 W SG z 







a 



ra 



2-v + vsg«g cos < 



G 



+ 8a 



Cj^aDaf) 2 
r3(af 3 -l)3 UsD - 



(16) 



Solving this simplified system produces results which agree quantitatively with the full 
system f TIU]) - f[T2^) for weak driving and qualitatively for all bias ranges. 

4. Static solutions and linear response 



We consider first the statics and small amplitude vibrations around equilibrium and 
derive a formula for the output power in this regime. To this end we expand ( |T0l) -( |T2l) 
to first order in v + /(v$ — Vq) and r/(«G,D — ^o), i- e -, we assume the oscillation amplitude 
to be small compared to the maximum amplitude allowed for a given static deflection 
and obtain 



7ru; 



{a G - x ) 2 
c° G a G v + v S G 

2o— r-r COS I 



r(a G - xqY 

Here xq is the deflection in the absence of drive, i.e., v + 



a 

Xq = — 
2 



c g «gWsg 



+ 



_(a G -a;o) 2 [a D - x ) 
and uu r is the renormalized frequency 



a 
2 



c g «gw| g 



+ 



(ctG-xar [a D -x ) 
We note that while the equation for xq can have two solutions only the solution with 
the smaller deflection is stable. The solution with large deflection is always unstable 
and leads to snap to contact. If surface forces are taken into account or if the drain 
electrode is placed outside the reach of tip bistable operation can be obtained as in [12J. 
From the above relations we find in the linear response regime a power transmission of 

,2 



\S- 



21 



4e 2 



_2. 2 .2 



G-D^ 2 



a; 2 ) 2 + ^ 2 uj : 



i 2 2 

+ UJ c GD 



+ 2aK G K B c GD - 



co 2 (u 2 -u; 2 ) 



with kq d 



(lu 2 — to 2 ) 2 + 7 2 Co> 2 
cg,d«g,d(ws - ^g,d)/(«g,d - ^o) 2 - 



;i7) 
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Figure 3. Typical shape of effective potential for the driven oscillator. Marked in the 
figure are the maximum x m above which the system is unstable towards snap-to-contact 
and the locus of the point where the second derivative changes sign, V^(x c ) = 0. 
Also shown are the typical oscillation regimes and their amplitudes. The qualititave 
difference between these regimes is best understood from looking at the phase (f> of the 
oscillator relative to the driving field (see figure O. 

5. Nonlinear response 

We now go beyond the linear response regime and consider the full solutions of ffT0]) - (fT2]) . 
For the general case these equations need to be solved numerically. In order to illustrate 
the typical resonant behaviour we will consider a specific system with a multi-walled 
carbon nanotube, D Q = 30 nm, D x = 20 nm extending a length L = 250 nm out from 
the support. The electrode dimensions are (see figured]) Zq = 100 nm, = 225 nm, 
Wd = Wt> = 50 nm, H = 25 nm and Hq = = 10 nm. Finite element modeling 
of this structure gives us the capacitances for an unbent configuration Cq = 5.4 aF, 
C D = 4.1 aF, Cqd — 6.1 aF and = 4 aF. For the mechanical properties of the 
nanotubes we have assumed an effective Young modulus of E = 1 TPa and a quality 
factor of Q = 200 (2?l EHl EHl [30l EH E21 E31 El] . 

For a qualitative understanding of the response it is useful to look at the non- 
averaged equations of motion. If we ignore excitations of other than the fundamental 
frequency we have the differential equation 

£ + 7^T + ^^eff(0 = <™+VSGCg(£) coswr 

with V e g = \ [£ 2 — o"^sd c d(0 — crt 'sG c G(0] • At a Source bias Vjg = 5 V and with 
V% = Vq = V this potential has the shape as shown in figure [3] and we can clearly 
distinguish a few different scenarios. For small excitations we expect to obtain the linear 
response solution stated above. As driving and amplitude increase we reach a point 
where deviations from the parabolic potential become manifest, leading to a hysteresis 
downward in the frequency plane. We expect this regime to have a similar frequency 
response as the Duffing oscillator. This regime is discussed below in subsection 5.1. As 
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Figure 4. Amplitude of oscillation in linear response (V+ = 0.03, 0.07 V) and 

hysteretic regime (V + — 0.09, 0.17 V). ( ) stable solutions, ( ) unstable 

solutions. 



one increases the driving further the oscillation amplitude is expected to increase. As 
it increases above the point where the curvature of the effective potential changes sign 
there will be a qualitative change in the response. This change is most markedly seen in 
the phase response (see figure E]) as "gap" around (f> = rr/2. We will discuss this regime 
in subsection 5.2. 



5.1. Onset of hysteresis 



For sufficiently small vibrations (fT0l) - ffl~2l) can be expanded in terms of Xq and r to obtain 
the frequency response equation 



(■yu) 2 r 2 + r 2 



(Co 2 - uo 2 ) - f3 2 r 2 



GCr 



«G 



2 „,2 



V , V, 



SO 



;is) 



where 



and 



2 



3<7(& , ac° D vi B + 2al 



4« D 



7 SD 



or 



uj 2 w 



a 



SG 



\(3 2 al. 



G 



This third order equation is the same as for the Duffing equation and is adequate for 
determining the onset of the hysteretic behavior. In the present case it performs less 
well to determine the frequency response for intermediate oscillation amplitudes in the 
hysteretic regime in which case one needs to solve either f|T0|) -f fT2|) or f|T4l) - f|T6l) . We 
illustrate this regime in figure H] where the amplitude of oscillation has been obtained 
from solving (|T0l-(|T2l) numerically The small inset of the figure illustrates the shape of 
the phase response across the resonance for this family of curves. The phase response 
for a larger set of voltages is shown in figure El 
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1.4 1.6 1.8 2 2.2 2.4 
Frequency [GHz] 

Figure 5. Relative phase of oscillations for varying AC-drive amplitude V+. Solid 
lines indicate stable solutions to (fTU)) - (fT2")) while dashed lines correspond to unstable 
solutions. For small drive amplitudes V+ < 0.05 V the system response agrees 
with linear response and only stable solutions exist. For intermediate values (V+ = 
0.07, 0.17 V) hysteresis in the frequency plane appears along high amplitude 
solutions (mainly unstable). Increasing drive amplitude further a gap in the phase 

response curve opens up around 4> = it/2 (curves V+ = 0.19, ...,0.27 V). ( ) stable 

solutions, ( ) unstable solutions. 

In figure [5] full lines correspond to stable solutions and dashed lines to unstable 
solutions. We note that there exists a set of low frequency solutions to the full equations 
Eqs. (jT0l) - (fl"2l) for the intermediate voltages (V + = 0.09, 0.17 V). These solutions 
cannot be described by the frequency response equation (fTBj) and are large amplitude 
solutions. The corresponding amplitudes and average displacements for these solutions 
are shown in figure [6) They are mostly unstable due to the fact that their peak oscillation 
amplitude exceeds x m . Note, however, (see figure [5] and figure [6]) that perturbation 
theory predicts that some of these solutions have stable (V + = 0.17 V) regions. These 
regions cannot be reached by sweeping the frequency down from a higher frequency but 
appears as a disconnected manifold of stable solutions. 

5.2. Large driving, large hysteresis 

As the driving increases even more the peak amplitude of the oscillator will eventually 
reach the point where the curvature of V e R changes from positive to negative. This is 
clearly seen in the phase response in figure [5] as a gap opens up around (p = ir/2. The 
corresponding large amplitude vibrations in this regime are rendered unstable when the 
oscillator amplitude reaches the point where it "rolls over the hill" (xq + r > x m ) and 
the system snaps into contact. The oscillation amplitude and average displacement for 
one such solution V + = 0.23 V are shown as the full lines in figure [7J 
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Figure 6. Large amplitude solutions in hysteretic regime (V+ = 0.09, 0.17 V). Most 
of the solutions have oscillation amplitudes too large to be stable. Note, however, that 
perturbation theory predicts a stable region for the bias V+ = 0.17 V (red curve). 
( ) stable solutions, ( ) unstable solutions. 

6. Validity of perturbation theory and output power 

We have also compared the analytical expressions with direct numerical integration 
of the differential equations. For small amplitude vibrations, including the primary 
hysteretic regime we find excellent agreement between numerical simulations and 
perturbation theory. Deviations from the predictions of perturbation theory are only 
seen for the large amplitude vibrations in the regime where there is a "phase gap". 
An example of such a comparison is shown in figure [7] where we clearly see the good 
agreement for small amplitude vibrations and the deviations at larger amplitudes. By 
looking at the detailed motion of the tube tip at large amplitudes one can see that the 
approximation by a pure harmonic motion is no longer a good approximation and higher 
harmonics have to be taken into account. 

Furthermore, perturbation theory predicts the existence of disconnected manifolds 
of stable large amplitude orbits that cannot be reached by sweeping down in frequency. 
We have not been able to detect any such orbits in the numerical simulations. This may 
be either due to the solutions being unphysical solutions to the perturbation theory 
equations or that appropriate initial conditions have not been used in the numerical 
simulations. 

The above only concerned matters related to the averaging method. Another issue 
is related to the assumption of large impedance mismatch (smallness of e). Typical 
transmission line impedances are of the order 50 Q. With capacitances in the aF range 
this means that we have e < 10 -6 which makes the approximation excellent. On the 
other hand this also means that transmitted power is very low. The obvious way to 
remedy this is to connect several devices in parallel, coupled to common gate and drain 
electrodes. For an array of N devices coupled this way one finds that Z should be 
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Figure 7. Oscillation amplitude r and average displacement Xq in the "phase gapped" 
oscillation regime. Solid lines represent result from perturbation theory for a bias of 
V+ = 0.23V while circles (O ) mark the result of numerical integration of the full system 
of ODE:s. For small amplitudes the agreement between numerical integration and 
perturbation theory is very good while perturbation starts to fail for larger amplitude 
due to the anharmonicity of the tip motion. 

replaced by an effective transmission line impedance of Zq G = NZ . Shown in figure M 
is the S^i-parameter for such an effective impedance of lOkQ calculated using ffTBT) for 
biases corresponding to the three different oscillation regimes. The arrows denotes the 
associated jumps in the hysteresis curve. 

Comparisons with numerical integration show that treating the system to lowest 
order in e as we have done here is a good approximation as long as |S 2 i| C 1. As 
| £21 1 ~ 1 the effect of dissipation in the tube-source contact and power delivered to the 
gate/drain result in a loaded Q-f actor exceeding the bare mechanical Q-f actor leading 
to a broadening of the resonances. 

7. Conclusions 

We have carried out an investigation of the RF-response in a three-terminal carbon 
nanotube resonator structure. By employing perturbation theory we have reduced the 
problem of determining and classifying the frequency response to that of solving a set of 
algebraic equations. We have found three distinct oscillatory regimes: linear response, 
hysteretic and a "phase gapped" regime which can lead to large hysteresis loops in 
the frequency response. Comparisons with direct numerical integration have shown 
that perturbation theory is qualitatively correct, and quantitatively accurate as long as 
the amplitude of oscillation is not too large. The perturbative treatment includes also 
terms related to parametric driving. We find no qualitative change in the behavior in 
the RF-response that arises from incorporating such terms. 
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Figure 8. |52i| 2 (= -Pin /-Pout) as a function of frequency for transmission line 
impedance of Z$ = lOkfl, or 200 devices in parallel connected to common 50 
transmission lines (see text). Note the disconnected stable manifold at V+ = 0.17V 
and /«1.8 GHz. 
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Appendix A. Derivation of ([8]) 

Here we outline the derivation of the lumped model for tube vibrations used in the 
paper. We begin by writing the solution to the PDE in the form 



oc 



u(x,t) = A /o (x) + ^2^ n (t)u n (x). 

n=0 

Here u n (x) are the eigenmodes with frequency Q n satisfying the homogeneous equation 
(-nl+(EI/pA)d*)u n {x) = while the offset A /o (x) satisfies (EI/pA)d*A h (x) = f {x) 
for some function fo(x) to be determined (typically this is the non- vanishing part of the 
time averaged force). Inserting the solution into the PDE for u(x, t) and projecting onto 
a normal mode u m (x) we find 

7m(*) +<£7m(*) = -L- 3 J o L dxul{x)Ux) + ^}^ t)]2 Jo dxul{%) 5 -^) 

We use here the normalization convention J L dxu^x^^x) = L 3 5 nm . We will typically 
work with a system where the driving only excites resonances of the fundamental mode 
uq(x). So we can quite safely make the approximation u(x,t) = A(x) + r yo(t)uo(x). The 
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differential equation then reads 

7o(t) + noTb(f) = y o d^(x)/ (x) + E d -^P" X rf -o ^ 

We will write this in terms of the displacement at the tip of the cantilever x^{t) = 
A(L)+7o(*)«o(£), 

x T (t) + figx T (t) = -^M [ L dxul(x)e(x) + Q 2 e(L) 

L 6 Jo 

+ 2^^ D Am Jo dXul{x) Su(xj • 

Here, the static bending has been expressed in terms of the deviation from the shape of 
the deflection in the lowest mode. A(x) = u (x)A(L)/u (L) + e(x). Assuming that the 
statically deformed shape closely resembles the shape of the resonance shape i.e., we set 
e(x) = 0, we arrive at 

Finally, using that u Q (x) = uq(x) one finds 



Mt) + nfar(t) = — E m(t)] 2 c t (x T ) 

2M eS i=G D 

(A.1) 



where M eS = ML 2 /u (L) 2 w M/5.684 and fi = 3.516L- 2 J EI / pA 
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